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1. Introduction 

One of the primary aims of quantum measurement theory is to under- 
stand the mechanism by which potential properties of quantum systems be- 
come actual. This is not an abstract or philosophical problem. Nowadays it is 
possible to carry out prolonged observations of individual quantum systems. 
These observations provide us with time series of data, and a complete theory 
must explain the mechanism by which these time series are being generated; 
it must be able to "simulate" the natural process of events generation. There 
are several methods of approaching this problem. John Bell Q for instance, 
sought a solution in hidden variable theories of Bohm and Vigier, his own 
idea of beables, and also in the spontaneous localization idea of Ghirardi, 
Rhimini and Weber 0. More recently, in a series of papers, two of us (Ph. 
B. and A.J.) |], |5| proposed a formalism that goes in a similar direction 
but avoids introducing other hidden variables beyond the wave function it- 
self. Our "Event Enhanced Quantum Theory" (in short: EEQT) describes 
a consistent mode of coupling between a quantum and a classical system, 
in which a classical system is one described by an Abelian algebra. It is an 
enhancement because it modifies quantum theory by adding the new term 
to the Liouville equation. This allows to unify the continuous evolution of 
a wave function with quantum jumps that accompany real world events. 
When the coupling constant is small, events are rare, and EEQT reduces to 
the orthodox quantum theory. 

We suggest that a measurement process is, by definition, a coupling of a 
quantum and a classical system, where transfer of information about quan- 
tum state to the classical recording device is mathematically modeled by a 
dynamical semigroup (i.e. semigroup of completely positive and trace pre- 
serving maps) of the total system. It is instructive to see that such a transfer 
of information cannot, indeed, be accomplished by a Hamiltonian or, more 
generally, by any automorphic evolutionQ To this end consider a system 
described by a von Neumann algebra A with center Z. Then Z describes 
the classical degrees freedom of the system. Let u be a state of A, and let 
u\z denote its restriction to Z. Let at be an automorphic time evolution of 
A, and denote uj 1 = a* (a;), where the dual evolution of states is given by 
a* (u) (A) = u>(a t (A)). Each a t is an automorphism of the algebra A, and so 

1 For a discussion of this fact in a broader context of algebraic theory of superselection 
sectors - cf. Landsman || Sec. 4. 4]. Cf. also the no-go result by Ozawa 
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it leaves its center invariant: at '■ Z — > Z. The crucial observation is that, 
because the evolution of states of Z is dual to the evolution of the observ- 
ables in Z, and we have a*(o;)|.g = o^lzi^lz), the restriction u^z depends 
only on uP\z- In other words the future state of the the classical subsystem 
depends only on the past state of that subsystem and - not on its extension 
to the total system. This shows that no information transfer from the total 
system to its classical subsystem is possible - unless we use more general, 
non-automorphic evolutions. The idea of describing a quantum measure- 
ment as a two-way coupling between quantum system and a classical system 
occurred before to several authors - we mention only the classical papers by 
Sudarshan || - but never within the completely positive semigroup approach. 

EEQT has several points of contact with other approaches. The mathe- 
matical model was a result of our studies of the papers of Jauch ||, [1(J, Hepp 
[|TJ, Piron Gisin |T§ |T§ and Araki Q , and also of the papers 



by Primas (cf. (18], [HJ). It was then found that our master equation describ- 
ing a coupled quantum-classical system is of the type already well known to 
statisticians. In his monographs [^, dealing with stochastic control and 
optimization M. H. A. Davis, having in mind mainly queuing and insurance 
models, described a special class of piecewise deterministic processes that 
was later found to fit perfectly the needs of quantum measurement theory, 
and that reproduced the master equation postulated originally by the two of 
us in IBfl. 



In [g2| it was shown that the special class of couplings between a classical 
and quantum system leads to a unique piecewise deterministic process with 
values on E'-the pure state space of the total system. That process consists 
of random jumps, accompanied by changes of a classical state, interspersed 
by random periods of Schrodinger-type deterministic evolution. The process, 
although mildly nonlinear in quantum wave function if), after averaging, re- 
covers the original linear master equation for statistical states. The action of 
the dynamical semigroup T t is given in terms of the process in the following 
way 

T t (P x ) = J P(t, x, dy)P y , 

where P(t,x,dy) is the transition probability function of the process and 
y — > P y is a tautological map, which assigns to every point y 6 E a one- 
dimensional projector P y . Let us discuss more precisely this connection be- 
tween a dynamical semigroup T t and the Markov- Feller process associated 
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with it. Suppose A4(E) is a Banach space of all complex, finite, Borel mea- 
sures on E with norm given by = \/i\(E). Let us define a map 



where At* is the predual space of the total algebra At- It is clear that ir is a 
linear, surjective and positive map with ||7r|| = 1. Therefore we can identify 
At* with the Banach quotient space M. (E) /ker7r. We say that a Markov 
process P(t, x, dy) with values in E is associated with T t iff U t = T t , where 
U t is the semigroup on A4(E) determined by P(t, x, dy), and U t denotes 
the quotient semigroup with respect to ker7r. The process is said to be 
Feller iff it preserves the space of all continuous and vanishing at infinity 
functions on E. It means that to find the associated process we have to 
extend the semigroup T t from M (E) /ker-zr to M. (E) in an invariant way. In 
general, such an extension may not exists or, if it exists, may not be unique. 
However, under mild assumptions, both the existence and the uniqueness 
were proved by analytical methods in ref. 22. From the point of view of the 
physical interpretation the uniqueness is of great importance since it leads to 
a unique description of the behavior of an individual quantum system under 
observation. 

The main objective of this paper is to provide a probabilistic construction 
of the process and discuss some of its properties and applications. We present 
a detailed construction since the theory of piecewise deterministic processes 
is not a part of the standard mathematical education. And, we believe, 
it is impossible to understand the essence of EEQT without having even a 
rough idea about this theory. The paper is organized as follows. In Section 

2 the formalism for classical-quantum interactions is presented. In Section 

3 the probabilistic construction of the PD process is described and some of 
its properties are analyzed. In Section 4 the classical part of the process is 
discussed. We also present an example of direct photodetection. Concluding 
remarks are given in Section 5. 

2. The Formalism 

We start by recalling the theorem by Christensen and Evans that de- 
scribes the most general form of a generator of a completely positive semi- 
group of transformations of an algebra with an non-trivial center. The the- 
orem generalizes the classical results of Gorini, Kossakowski and Sudarshan 
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ZJJ and of Lindblad [51 to the case of arbitrary C*-algebra, and it states 
that essentially the Lindblad form of the generator holds also for this more 
general case. We quote the theorem for the convenience of the reader [EH]: 



Theorem 2.1.(Christensen - Evans) Let a t = exp(Lt) be a norm- continuous 
semigroup of CP maps of a C*- algebra of operators A C B{7i). Then there 
exists a CP map <fi of A into the ultraweak closure A and an operator K £ A 
such that the generator L is of the form: 

L(A) = <f>(A) + K*A + AK 

Let us apply this theorem to the case of A being a von Neumann algebra, 
and the maps a t being normal. Then can be also taken normal. We also 
have A = A, so that K £ A. Let us assume that a t (I) = I or, equivalently, 
that L(I) = 0. It is convenient to introduce H = i(K — K*)/2 £ A, then 
from L(I) = we get K + K* = -</>(/), and so K = -iH - 0(1) /2. Therefore 
we have 

L(A)=i[H,A] + <l>(A)-{<l>(I),A}/2 

where { , } denotes anticommutator. 

We now apply the above formalism to the hybrid system which is a di- 
rect product of the classical and quantum mechanical one. The physical idea 
behind such a model is that a quantum measurement is to be defined as a 
particular coupling between a quantum and a classical system. We continu- 
ously observe the classical system, notice changes of its pure states (we call 
these changes "events") and from these we deduce properties of the coupled 
quantum system. Details can be found in refs. 4 and 5. One can think 
of events as 'clicks' of a particle counter, sudden changes of the pointer ve- 
locity, changing readings on an apparatus LCD display. The concept of an 
event is of course an idealization - like all concepts in a physical theory. Let 
us consider the simplest situation corresponding to a finite set of possible 
events. The space of pure states of our classical system C, denoted by S c , 
has m states, labeled by a — 1, . . . , m. Statistical states of C are probability 
measures on S c - in our case just sequences p a > 0, J2 a Pa = 1- 

The algebra of observables of C is the algebra A c of complex functions 
on S c - in our case just sequences f a , a = 1, . . . , m of complex numbers. We 
use Hilbert space language even for the description of the classical system. 
Thus we introduce an m-dimensional Hilbert space 7i c with a fixed basis, 
and we realize A c as the algebra of diagonal matrices F = diag(/i, . . . , f m ). 
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Statistical states of C are then diagonal density matrices diag(pi, . . . ,p m ), 
and pure states of C are vectors of the fixed basis of ri c . Events are ordered 
pairs of pure states a — > (3, a ^ f3. Each event can thus be represented by an 
m x m matrix with 1 at the (oc,f3) entry, zero otherwise. There are m 2 — m 
possible events. Let us point out that important here is the discreteness of the 
classical system not its finiteness. We can easily generalize the above to the 
case when the classical points form, for example, the set of natural numbers. 
Then the classical algebra becomes l°° (uniformly bounded sequences) while 
statistical states are positive elements from I 1 (summable sequences). 
We now come to the quantum system. Let Q be the quantum system whose 
bounded observables are from the algebra A q of bounded operators on a 
Hilbert space ri q . In this paper we will assume ri q to be finite dimensional 
. Pure states of Q are unit vectors in H q ; proportional vectors describe the 
same quantum state. They form a complex projective space CP(ri q ) over 
ri q . Statistical states of Q are given by non-negative density matrices p, 
with Tr (p) = 1. 

Let us now consider the total system T = Q x C . For the algebra At of 
observables of T we take the tensor product of algebras of observables of Q 
and C: At = A q ® A c . It acts on the tensor product Tt q ®ri c = ©™ =1 7i Q , 
where 7i a ~ 7i q . Thus At can be thought of as algebra of diagonal m x m 
matrices A = (a aj 3), whose entries are quantum operators: a aa e A q , a a @ = 
for a ^ f3. Statistical states of Q x C are given bymxm diagonal matrices 
p = diag(pi, . . . , p m ) whose entries are positive operators on H q , with the 
normalization Tr (p) = X^Tr (p Q ) = 1. Duality between observables and 
states is provided by the expectation value < A > p = J2 a Tr {A a p a ). 

We will now generalize slightly our framework. Indeed, there is no need 
for the quantum Hilbert spaces TC a , corresponding to different states of the 
classical system, to coincide. We will allow them to be different in the rest 
of this paper. Intuitively such a generalization corresponds to the idea that 
a phase transition can accompany the event. We denote n a = dim(ri a ). 

We consider now dynamics. It is normal in quantum theory that classical 
parameters enter quantum Hamiltonian. Thus we assume that quantum 
dynamics, when no information is transferred from Q to C, is described by 
Hamiltonians H a : ri a — ► H a , that may depend on the actual state of 
C (as indicated by the index a). We will use matrix notation and write 
H = diag(i? Q ). Now take the classical system. It is discrete here. Thus it 
can not have continuous time dynamics of its own. 
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The coupling of Q to C is specified by a matrix V = (gap), where g a p 
are linear operators: g a p : Tip — > 7i a . We assume g aa = 0. This condi- 
tion expresses the simple fact: we do not need dissipation without receiving 
information i.e without an event. To transfer information from Q to C we 
need a non-Hamiltonian term which provides a completely positive (CP) 
coupling. As in 0, [5| we consider couplings for which the evolution equation 
for observables and for states is given by the Lindblad form: 

A a = i[H a , A a ] + 9p a A p9pa ~ ^{A Q , A a }, 

or equivalently: 

p a = -i[H a ,p a ] +Y.9apppg af3 ~l-{A a ,pa}, 

13 l 

where 

Aa = 9%9pa 

P 

The above equations describe statistical behavior of ensembles. Individual 
sample histories are described by the following algorithm: 
Suppose that at time to the system is described by a normalized quantum state 
vector ip and a classical state a. Then choose a uniform random number 
p G [0, 1], and proceed with the continuous time evolution by solving the 
modified Schrddinger equation 

tpt = {-iH a - ]^A a )i)t 

with the initial wave function ipo until t = t\, where t\ is determined by 

h 

J (ip t , A a ijjt)dt = p 

to 

Then jump. When jumping, change a — > (3 with probability 

p a ^ p = \\gpaipti\\ 2 /(Ai, A a ip tl ) 

and change 

tpu = gpa^tJ\\gpaAi\\- 
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Repeat the steps replacing to, ipo, a with t\, ipi and (5. 

This leads to a stochastic process, in which the randomness appears as 
point events i.e. there is a sequence of random occurrences at random times 
T\ < T 2 < but there is no additional component of uncertainty between 
these times. It consists of a mixture of deterministic motion and random 
jumps. A class of such processes is called piecewise deterministic processes 
(PDP) [ffifl . The motion between jumps is determined by a complete vector 
field X on the pure state space E of the total system. The jump mechanism 
is determined by two further components: a jump rate A and a transition 
kernel Q. The vector field X generates a flow <p(t, x) in E, which is given 
by <p(t, x) = J x (t), where ■y x (t) is the integral curve of X starting at point 
x G E. The jump rate is a measurable function A : E — > R + such that for 
any x G E the mapping t — > A o (f)(t, x) is integrable at least near t = 0. 
The set of those x G E for which X(x) = we denote by E . The transition 
kernel Q : 13(E) x E — > [0, 1] satisfies the following conditions: 

a) Q(E, x) = 1 Vx G E, 

b) Q({x}, x) = if x e E\E and x) = 1 for x G £ , 

c) VT G B(E) the map x — > Q(T, x) is measurable. 

Here 13(E) denotes the Borel cr-algebra on E. In our case E = \JCP a , a = 
1, 2, ...,m and we have the following formulas for X, A and Q: 

X/(^, a) = — /(- — \ a) t=0 

A( , 0, a) = < if), A a ip > 

Q(d0, /3; V, a) = tt; rS(<f) - t, -rMcf) 

Kv>> a ) \\9fiaVW 

if (if), a) G E \ E and S((f))d(f> denotes the Dirac measure. 
The triple (X, A, Q) is called local characteristic of the process. Its infinites- 
imal generator is given by 

Cf(x) = Xf(x) + X(x) J [f(y) - f(x)]Q(dy, x) 

E 

and produces sample paths exactly such as described by the above algorithm. 
It is worth noting that there are no correlations between jump times T{, i G N. 
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However, as we will see in the next section, the survival function of random 
variable T\ is exponential, given by 

t 

F t = exp(- J X(x s )ds) 
o 

where x s denotes the piecewise deterministic process. 



3. The PD process 

In this section we present the detailed construction of the process intro- 
duced in Section 2 and investigate some of its properties. General references 
on stochastic processes are [E7L 123]. Probabilistic concepts can be found in 



2S, 30 



At first we construct a probabilistic space (O, A) (compare [|3T| for a 
similar construction for Markov decision processes). Let Q be a set of all 
sequences (t , x ; £i> Xi\ ■ ■ ■), which are finite or infinite, and such that t = 0, 
t n < t n +i, t n G R+ = [0, oo], x n G E for all n G N U {0}. If a sequence is 
finite i.e. uj = (t , xq; . . . , t n , x n ) then we put 

t n +l — t n +2 — ■ ■ ■ — OO, X n +i = £ n +2 — • . . = X n 

oo 

It follows that Q can be embedded into an infinite product space n 

n=0 

where fio = {0} x E and Q n = R + x E 1 . On each Q n we have a natural 
a-algebra ^4 n given by £>(R+) ® B(E). We define a a-algebra ^4 on Q as 
(®^L X)|n. 

Now let us construct a family of probabilistic measures P x on (fi, ^4) with 
respect to an initial state x G E. They will be determined by the determin- 
istic drift 0, the jump rate A and the transition kernel Q. Because we want 



to use the Ionescu Tulcea theorem [29| we have to define transition kernels 



between (fl n , A n ) and (fi n+1 , A n+ i). We do it step by step. 

On tt we take the Dirac measure P = 8 X . Let A(t, x) := Jq A(0(s, x))d& 

and let us define 

F x {ti) = 1 -exp(-A(t 1 , x)) 
K x (ti, dxij = Q(dxi, <j)(ti, x)) 
As the transition kernel between (Q ,Ao) and (Qi,Ai) we take 

P^(x, B x x ri) = If K x (t u dx 1 )dF x {ti) 
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for any B\ G £>(R+) and any Y\ G B{E). In the second step we define 

f a \ - J °> if *i > t2 

^WO^J - li-expC-A^-^xO), if*!<* 2 
K(t uxi )(t2, dx 2 ) = Q(dx 2 , (f>(t 2 -h, xi)) 

and put 

P 2 1 (*i> zi; B 2 xT 2 ) = J J K {tljXl) (t 2 , dx 2 )dF {tuXl) {t 2 ) 

It is clear that P 2 X is a transition kernel between .Ai) and (fi 2 , A 2 ). In the 

similar way we construct higher kernels P" +1 . By Ionescu Tulcea theorem 

00 

there is a unique probabilistic measure P x on ( ^n, ®^L .4. n ) sucn that 

n=0 

for every measurable rectangle A = A x A\ x . . . x A n x fi n+ i x . . . the 
following identity 

P X [A] = 5 X (A ) J Pq(x; dti, dxt) ■■■ j i^ _1 (t„_i, x n _i; dt n , dx n ) 

is satisfied. It is clear from the above formula that P x is concentrated on 
il x = {lu £ VI : x = x}, x G E. Moreover P x is measurable with respect to 
x. 

To investigate properties of the above measure let us define a sequence of 
measurable random variables 

T n : Q. x — > R+ T n (co) = t n , X n : Q x — > E X n (u) = x n 

The distributions of T and X are Dirac measures concentrated in {0} and 
{x} respectively. The distribution dF Tl of Ti is given by 

P x [Ti < t] = 1 -exp(-A(t, x)) 

and the conditional expectation of X 1 given T x equals to 

S x [l {XlC r}|7i] = Q(T, 0(Tx, x)) 

Here 1{.} denotes an indicator function of a given set. The above equation 
can be also written as 

dF Xl \ Tl (y\t) = Q(dy, (f)(t, x)), 
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where the left hand side is the conditional distribution of X\. For arbitrary 
nGNwe have the following formulas: 

w n IT Y 1 - / if t < T n 

ArLl{T„ +1 <t}|Jn, ^nj - \ 1 _ ex p(-A(t - T n , X n )) if t > T n 

E x [l{x n+1 dT}\X n , T n+1 ] = Q(T, 0(T n+ i, X n )) 

It follows that P x [Ti = 0] = so T\ > a.s. Because, after a jump, process 
starts again so T n < T n+ \ a.s. for every n. This fact can be also derived from 
the following equality: 

P x [T n+1 -T n >s] = £ x [exp(-A(s, X n ))\ 

It means that a set of paths with two or more simultaneous jumps has zero 
probability. Moreover, because Q({x}, x) — for every x G E \ E so with 
probability one the process can not jump to the state it is deterministically 
approaching. There are no jumps from the set E at all. 
Let us calculate some physically interesting probabilities. For example the 
probability that there is no jump up to time t equals to 



z 

P^ > t] = exp(- J A(0(s, x))ds) 



Because 

P X [T 2 > t] = P x [7\ >t] + £7 x [l {Tl < t} exp(-A(* -T u X ± )] 
and, on the other hand, 

P X [T 2 >t] = P x [T 2 >tAT 1 <t] + P X [T 2 > t A Ti > t] 

so the probability that exactly one jump happens up to time t is given by 

t 

P X [T 2 >tAT 1 <t] = J l {u < t} dF Tl {u) J E exp(-A{t - u, y))dF Xl \ Tl (y\u) = 

o 

/ / \((f)(u, x)) exp(— A(u, x)) exp(— A(t — u, y))Q(dy, <f>(u, x))du 

JO J E 
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Now let us define a random variable = lim^oo T n . For every t < we 
construct the process x t by putting 

x t (w) = <f>(t - T k (u), X k (u)) if T k {u) < t < T k+1 (cu) 

In general we can have the process with the lifetime. We show that in our 
case, due to the boundedness of the jumping rate, = oo a.s. Let C = 
sup xeE X(x). Then for every t > 

sup(l - exp(-A(i, x)) < 1 - e~ ct 
xeE 

Let us fix t and denote C\ — 1 — e~ ct , which is strictly less than I. Then 
P x [T n+1 < t] = E[l {Tn < t} (\ - exp(-A(i - T n , X n ))] < C x P x \T n <t]< C? +1 
by induction. It implies that 

oo 

P x [f] { T n < t}} = lim P x [T n < t] = 

n=0 

It follows that xt is defined for all t G R+ and is a cadlag process i.e. pos- 
sesses right continuous with left limits paths. 

To end the construction of ingredients needed for a Markov process let us 
introduce a natural filtration on Q x given by = <r{x s , s < t} and take 
= VtFt- Let T t and denote the P^-completion of and re- 
spectively. Because, after a jump, the process evolves deterministically, so 
the filtration (J r t)o<t<oo is right continuous. Hence the filtered probability 
space (Q, Too-, P x , Ft) satisfies the usual hypothesis for every x G E and 
x t is an adapted and cadlag process. Because the distribution of Ti de- 
pends only on the current state x t and, after a jump, process starts again so 
(Q, Too, P x , J-f, x t ) is a strong Markov process with infinite lifetime. 

Next we show another important property of the process x t , namely the 
quasi-left-continuity. Let us define a random set A = {(t,cu) : x t - ^ x t }, 
where x t - is the left limit of x t . Then 

v(u\ dt, dx) = 5^1a(s, w ) 5 (s,x s H)(^, dx )i 
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where 8r s ,x) is the Dirac measure on R + x E concentrated in (s, x), is an 
integer- valued random measure. It leads to a simple point process N t given 
by 

oo 

Nt = i/([0, t] x E) = £ l {Tn < t} 

71=1 

Because = oo a.s. so iVj is a.s. finite valued. It is also integrable because 

oo oo 
E x [N t ] = $>P s [T n <t]< J2 nC l < 00 



n=l n=l 



Moreover it was shown in p6| that the compensator of N t is equal to J * A(x s )ds 
and so M t := Nt — J A(x s )ds is an (P^, ^) -martingale. Using this fact it 
can be calculated that the dual predictable projection of v is given by 

u p (uj;dt, dx) = Q(dx, x. t (uj))\(xt(u))dt 



Thus v p (uj; {£}, E) — and so x t is quasi-left-continuous P§| . Thus we 
proved that x t is a Hunt process. Moreover x t is a Feller process i.e. the 
transition kernel of x t generates a strongly continuous semigroup of contrac- 
tions on the space of all continuous functions on E, see |3"3" |. 

4. Stochastic representation of the classical system 

In this section we discuss some properties of the stochastic process asso- 
ciated with the measuring apparatus. Let C be a state space of the classical 
system i.e. C = {1, 2, ...,m}. Let us define a {0, l}-valued process pf by 

P?(w) = <^ (xtH) , 

where 5p is the Kronecker delta and 7r denotes the canonical projection 71 : 
E —>■ C. By P(t, x, T), x E E and r C E we denote the transition kernel 
of the process x 4 . It was shown in |22|] that P(t, x, T) is associated with the 
dynamical semigroup T t and so 



T t (P x ) = J P y P(t, x, dy) 



Here P y is the one-dimensional projector corresponding to y G E i.e. P y = 
\y >< y\. We show that the average of pf gives the probability of finding 
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the total system at time t in a classical state a. Let = E[p^}. Then 
P? = J K(y)P(t, x, dy) = J Tr(P y )P(t, x, dy) = 

E CP a 



Tr( J P y P(t, x, dy)) = Tr(T t (P x 



CP a 

Now we derive a differential equation for p". Let us start with the following 
example. 

Example 1. Let C = {1, 2}. Then a change of the process p", a — 1, 2, is 
given by 

dp] = -p]-dN t + p 2 t -dN t , p\ = 1 

dpi = -p 2 t -dN t + p]-dN t , pi = 0, 

where N t is the counting process introduced in the previous section. Solving 
the above equations we get 

! + (-!)*' 1 ~ (-1)*' 

Pt ~ 2 , Pt - 2 

Because M t = N t — J" * \(x s )ds is a martingale so we get the following 
equations for averages p": 

dp] = E[(-pl+pl)\(x t )]dt 

df t = E[(-p 2 t +pl)X^ t )]dt 
When the intensity is a constant function equal to A they reduce to 

dp l t = \{-p\+p 2 t )dt, dp 2 t = X{-pl+p])dt 

with solutions given by 

_! 1 + e' 2Xt _ 2 1 - e~ 2Xt 

Pt ~ 2 , Pt ~ 2 

Proposition 4.1. 

dp? = -E[\ a (x t )]dt + Yl E[\\g a ^ t f\dt 
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Proof: Because pf = Tr(p a ), p a = T t (P x ) a so 



= Tr(p Q ) = Tr(-i[H a , p a \ - \{K a , p a } + 



dt ^ V -L"«, r«j 2 

S 9a/3p/3g*a(3) = Tr(A a )p a ) + ]T Ti{g* ap g a pp p ) 

f3+a f3+a 

On the other hand 

£[A a (x t )] = / < 2/|A a |2/ > P(t, dy) = Tr^T^),) 

CP a 

E[\\g a px t \\ 2 ] = J < y\g* af3 g af3 \y > P(t, x, dy) = Tr(g* af3 g a pT t (P x ) p ) 

CP,3 

so the assertion follows. □ 

The advantage of this stochastic representation of Tr(T t (P x ) a ) is that we can 
predict the future of the classical system if we know its past. Let us point 
out that the classical component of x 4 usually is not a Markov process. 

Let us assume that we start at t — with a quantum state x G CP a o, 
and up to the present we have observed the following classical trajectory 

(t = 0, a ), (h, ai), (t k < t, a k ) 

Then the probability p a that the next jump will go to a can be obtained as 
follows. Let us calculate 

X) fi'a fe a fc _ 1 0(4 — tk-1, Xk-l) 

•^1 II / / , \ || 5 "^h 



\9ctiao <f>(tl, x )\\ Wgakak-i&itk — tk-i, X k -l)\\ 

and x t = 4>(t — tk, Xk) for t > tk- Then 

oo 

Pa = E[\\g aak (x Tk+1 )\\ 2 } = J dF Tk+llTkjXk (t\t k , x k )\\g aak (x t )\\ 2 = 



oo ( 

J dtexp(- J \(x s )ds)\\g aak (x t )\\ 2 

These probabilities can be also used to determine an initial quantum state. 
Let us assume that we start at t — with a classical index a and with one of 
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the following pure quantum states x % G CP a o, i = 1, 2, n. The probability 
that the first jump will change the classical index onto a is given by 

Pa = £[||^ Q0 (x Tl )|| 2 ], T = 0, X = X* 

In the similar way we calculate probabilities p a2a i that the first jump will 
go to cti and the second one to a 2 and so on. Taking appropriate gp a we 
can make these probabilities significantly different for each initial quantum 
state x % and thus conclude which one is the most probable by observing the 
classical trajectory 

Example 2. Let us consider the fluorescent photons emitted by a single, 
two-level atom that is coherently driven by an external electromagnetic field. 
It is known that the quantum system evolves from the ground state in a 
dissipative way. When a photoelectric count is recorded by a photoelectric 
detector (we assume the detector efficiency to be equal to one), the atom 
returns to the ground state with the emission of one photon. Thus, after the 
emission of each photon, the atom starts its evolution from the same state. 
We describe this situation using the probabilistic framework introduced in 
the previous sections. 

The quantum system as a two-state system is represented by 2 x 2 complex 
matrices. The classical system, which counts emitted photons we describe 
by an infinite sequence of numbers n = 0, 1, 2, ... Hence the state space of 
the total system is equal to 



E = (J CP 2 



„P 2 

n=0 

The time evolution of the quantum system is described (for every classical 
index n) by the modified Schrodinger equation 

tpt = -iHipt = {-iH - -A)ip t , 
where A = / yA*A and 

The coupling operators g nm are given by g n +i, n = A and g mn = if m ^ 
n+ 1. A solution for ip t can be written as ipt = U(t)i/;o, where pi] 

•,/> „ rf /„ / cos/it — -r-sinut ZTp-sinut \ 

U(t) = e~ itH = e" 7t/4 n , 2/i , 7 • , 

y sin fit cos fit + ^- sin /it J 
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Here 7 is the relaxation rate, Q is Rabi frequency and ji = ^fl 2 — {j/2) 2 . 
The ground state ipo is given by 



^0 = 

The deterministic flow is defined by 

(f>(t, (ip , n)) 




1 



\\u(t)M 



for every n G N U {0}. The jump rate A is given by A((^, n)) = < ip, Aip > 
for all n and the transition kernel 

Q(d(f>, to; ip, n) = - ipo)d<f) 

Because of the uniqueness of a jump after the emission of a photon the 
classical component of the piecewise deterministic process of the total system 
is also a Markov process. Let us derive the distribution of the waiting time 
between jumps. In this case it is exactly the distribution of the random 
variable 7\. Thus 

t 

F(t) = 1 - exp[-A(f, (V> , 0)) = 1 - exp[-/ A(0(s, 0))ds] 



= 1 - exp[- / < gO* , ;A-A ^ > *] 

Because 

^-\\U(s)^ \\ 2 = < -iHU(s)i/j , U(s)^ > + < U(s)iJ>o, -iHU(s)^ > 
as 

= < U(s)^j , t(H* - H)U(s)^ > = - < U(s)^ , 1 A*AU(s)^ > 
so we obtain that 

r d 

F(t) = 1 - eMj ( T J°g\Ms)4> \\ 2 )ds] = 1 - \\U(s)4> \\ 2 
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Its density is given by 

f(t) = 7\\AU(t)i/>o\\ 2 = 7^sin 2 (/it)exp(- 7 t/2) 
It is exactly the waiting time density obtained in [fj5|| . 

Now let us consider the time evolution of the averages of the classical com- 
ponents of the process p n (t). In the present context they have a simple inter- 
pretation: p n (t) is the probability that n photoelectric counts are recorded 
in the time interval [0, t}. Hence p n (t) = P[T n < t A T n+i > t] and so 

p (t) = \\u(t)M 2 

t 

Pi(t) = [p (s)f(t-s)ds = (p */)(*) 



t t—S2 

p 2 {t) = J ds 2 p (s 2 ) J /(si)/(t - s 2 - sjdsi = (p *f*f)(t) 



and so on. In the above we extended functions p (t) and f(t) to the whole 
real line (— oo, oo) by putting value zero for negative arguments. The sign * 
denotes the convolution. Taking the Laplace transform of p n (t) with respect 
to variable t we obtain that 

oo 

p n (X) = fe- xt p n (t)dt = R,(A)[/(A)r 



which coincides with the formula given in (see also |[34}| ). 

The above simple example demonstrates how the piecewise deterministic 
dynamics works in practice. It also shows that EEQT reproduces known re- 
sults in certain domains. But it also predicts more in other domains, where 
the orthodox quantum theory is silent. For example the EEQT algorithm has 
been applied to tunneling time problems f36| . In particular, traversal and 



reflection times of electrons through a one-dimensional barrier have been 



calculated in pf\ . There EEQT gives predictions which can be tested and 
compared with those stemming from other approaches. In this respect, the 
orthodox quantum theory gives no prediction at all. It is also possible using 
the framework of EEQT to consider a simultaneous measurement of noncom- 
muting observables |38| . In the simplest case of a simultaneous measurement 
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of several spin projections for a spin 1/2 quantum particle, the piecewise 
deterministic process turned out to be a non-linear version of the Barnsley's 
iterated function system and led to the chaotic behavior and fractal structure 
on the space of pure states . 



5. Concluding remarks 

The crucial concept underlining our approach to quantum measurement 
is that of a classical and irreversible event. This is taken into account by 
including from the beginning classical degrees of freedom. From the struc- 
tural point of view such a coupling (EEQT) consists of the following essential 
ingredients: 

- tensoring of a non-commutative quantum algebra of observables with a clas- 
sical commutative algebra (or, more generally, taking the classical Abelian 
algebra as the center of the total algebra of observables), 

-replacing Schrodinger unitary dynamics by a completely positive semigroup 
describing the time evolution of ensembles, 

- interpreting the continuous time evolution of statistical states in terms of 
a piecewise deterministic process with values in the pure state space of the 
total system, 

- applying the uniqueness theorem for deducing the piecewise deterministic 
algorithm generating sample path of an individual system. 

This (EEQT) gives a minimal extension of the quantum theory which ensures 
the flow of information from the quantum system to the classical variables. 
Further, EEQT provides a way to calculate numbers needed in real experi- 
ments and also allows for natural mathematical modeling of feedback during 
experiments with quantum systems. 

This coupling with classical variables does not mean we are taking a step 
backward into classical mechanics. We are only saying that not all is quan- 
tum and there are elements of nature that are not, and clearly cannot be, 
described by a quantum wave function. Even if this is viewed as an assump- 
tion, it is firmly confirmed by experiments which show that we are living 
in the world of facts, not in the world of possibilities. Thus, for this as- 
pect, which is clearly not reducible to quantum degrees of freedom, we have 
adopted the term "classical variables." However, this does not imply that 
we intend to impose any strict restrictions on their nature; they may prove 
to be related, for example, to gravity or even to consciousness. If there is 
no interaction between these classical variables and quantum degrees of free- 
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dom, both systems may prove to evolve separately according to their own 
equation of motion. It may be that when they interact, a new dissipative 
operation appears which results in the irreversibility of the evolution and 
leads to collapses in measurement situations. It is worth noting that there is 
a fundamental difference between the classical variables in EEQT and addi- 
tional parameters in hidden variable theories. Hidden variable theories use 
microscopic variables that are hidden indeed from our observations. EEQT 
deals with classical variables that can be observed. In fact, it states that 
there are no other variables that can be directly observed. They are a direct 
counterpart of physics on the other side of the Heisenberg-von Neumann cut. 
Further, in hidden variable theories there is no back action of these variables 
on the wave function. In EEQT there is such an action. 

Let us also discuss two possible generalizations of the above framework. 
The first one concerns the dimension of the quantum Hilbert space. Here, for 
the sake of simplicity, we have used only finite dimensional Hilbert spaces, but 
from the construction it is quite clear that this assumption is not essential. 
We can admit infinite dimensional Hilbert spaces, but this requires additional 
work so that the formulas are well defined. For example, the infinite series 
of operators must be convergent. Further, we can allow the Hamiltonian 
operator to be unbounded. Also the existence of the deterministic flow can 
be established since CP(7i q ) is an infinite dimensional Hilbert manifold; that 
is, it can be covered by a family of open sets each of which is homeomorphic 
to an open ball in a Hilbert space. 

The second generalization is connected to the discreteness of the classical 
system. Although the basic applications of EEQT concerns measurement 
processes, it can also encompass a non-trivial interaction between a quantum 
system and a classical continuous one. The most transparent example is the 
SQUID-tank model, which consists of an electric oscillatory circuit coupled 
via a mutual inductance to a superconducting ring. In that system the 
oscillatory circuit acts as an external flux source for the SQUID ring, which 
induces a screening current in the ring. This screening current is coupled 
back to the classical circuit due to the mutual inductance. It results in the 
modification of the differential equation for the damped classical harmonic 
oscillator by the expectation value of the superconducting screening current 
operator. Now the classical phase space is a symplectic manifold R 2 and the 
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equation of motion is given by 

C0 + ^ + j<f> = I{t) 

where <fi is the classical flux variable. The quantum object evolves according 
to Hamiltonian 

g 2 ($-$,) 2 $ 

^ = 2C + 2A ^os( 27 r-) 

where Q = — i-^ is the momentum operator, $ is the position operator, 



$o — 2e anc ^ ^ ^ s ^ ne ex ternal magnetic flux. Suppose p t is an evolving 
statistical operator of the total system and define the so-called collective 
classical variable 

x t — I x(Trp t (x, p))dxdp 
Then its evolution equation takes form (we omit the physical constants): 



% + % + (1 + \\f\\h)x t = I(t) + \\f\\ 2 L2 <I S > 



p 



where I s is the screening current operator, / is a function monitoring the 
strength of the coupling and p is the reduced density matrix obtained by 
tracing over classical variables. Here < I s >p denotes the expectation value 
of I s in state p. For more details see ref. 40. It is worth emphasizing that 
the above modification of the evolution of the collective classical variable, 
which was postulated by experimental physicists (see eq. (12) in ref. 41), 
has been derived in a rigorous mathematical way within the framework of 
EEQT. Finally, using EEQT, a possible influence of the quantum matter 
on the classical gravitational field has been demonstrated. It was achieved 
by the change of a dynamical path of the classical particle moving freely 
along a geodesic curve when interacting with the quantum system. In av- 
erage the classical evolution equation is perturbed by the expectation value 
of the quantum position operator j|2| . The probabilistic description of the 



dynamics in this case was presented in ref. 43 and the back action of the clas- 
sical variables on the quantum system, resulting in its non-trivial asymptotic 
behavior, was discussed in ref 44. 

To sum up: EEQT, although it is a phenomenological model, is very 
promising and can be successfully applied to a large class of physical phe- 
nomena. 
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